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A classical continuum model of an effectively one-dimensional ferromagnet with exchange and anisotropics of hard and easy- 
axis type is considered. If an external field is applied along the easy axis, the lowest lying topological excitations are shown 
to be untwisted or twisted pairs of vr-domain walls. The fluctuations around these structures are investigated. It is shown 
that the fluctuations around the twisted and untwisted domain wall pair are governed by the same set of operators. The 
untwisted domain wall pair has exactly one unstable mode and thus represents a critical nucleus for magnetization reversal 
in effectively one-dimensional systems. The twisted domain wall pair is stable for small external flelds but becomes unstable 
for large magnetic flelds. The former effect is related to thermally induced coercivity reduction in elongated particles while 
the latter effect is related to "chopping" of twisted Bloch wall pairs in thin films. In view of a statistical mechanical theory 
of magnetization reversal which will be presented in a separate article, the scattering phase shifts of spin waves around these 
structures are calculated. The applicability of the present theory to magnetic thin films is discussed. Finally, it is noted that 
the static properties of the present model are equivalent to those of a nonlinear cr-model with anisotropics and an external field. 

PACS numbers: 03.50.Kk, 75.10.Hk, 75.60.Ch, 75.70.Kw 



I. INTRODUCTION 



Macroscopic ferromagnetic sampleJ3'i consist of many 
domains in which the magnetization is uniform and 
directed along one of the minima of the crystalline 
anisotropy. The phase boundaries between such regions 
are formed by domain walls (Bloch walls) in which the 
magnetization vector rotates continuously between differ- 
ent anisotropy minima. The formation of these domains 
is due to the long-range magnetostatic forces which tend 
to avoid the formation of magnetostatic charges at the 
sample surface. However, domain walls have locally pla- 
nar symmetry and can therefore locally be described by 
an effectively one-dimensional modeltrH. 

A one-dimensional description is also adequate for 
elongated samples of mesoscopic size if the lateral sample 
extension is less than a domain wall width. Such particles 
are wideli^ used in magnetic recording media, e.g. Cr02 
-particlesH are amost perfect needles with aspect ratios 
of up to 20:1. For this reason and in view of tremendous 
recent progress in sample preparation on the nanometer 
scale, it is therefore of particular importance to study 
the model of an effectively one dimensional ferromagnet 
in detail. 

In the following we shall focus on a description of the 
magnetization within a classical field theory. Such a for- 
mulation also provides the starting point for a quantum 
mechanical theory in the semiclassical limitu. The mag- 
netization is treated as a classical vector of constant mag- 
nitude and adjacent moments interact via exchange thus 
giving rise to a "stiffness" of the spin chain. The present 
model contains single ion anisotropics of hard- and easy- 
axis type which may have demagnetizing or crystalline 
origin. In addition it includes an external field directed 



along the easy-axis. Without an external field, this model 
is also known as a "biaxial ferromagnet" . 

The present model has also been used to describe 
weaklAL coupled one-dimensional (ID) ferromagnetic 
chainsOpi In effectively ID antiferromagnets such as 
TMMCeI, it emerges as an effective model for the sublat- 
tice magnetization. The dynamic version of this model 
without external field and damping has been shown to 
be integrableCI and rjeveals a surprisingly rich palette of 
soliton and breatherQ solutions, the solitons playing the 
role of domain walls. 

The simplest static, topological excitation in a biax- 
ial ferromagneLin the absence of an external field is the 
TT-Bloch walloEj (see FigJ|) which constitutes the tran- 
sition region between two equivalent anisotropy minima. 
Its stability has been investigated by Winterliil who ex- 
plicitly derived spin wave excitations. He showed that 
within the ID system, Bloch walls are stable save for the 
zero energy mode which describes a rigid translation of 
the domain wall. Later Janakt3 quantized the spin wave 
excitations around a pinned domain wall and included 
demagnetizing effects of spin waves running parallel to 
the domain wall. Hornreich and ThomasE3 considered 
a biaxial ferromagnet with an external field perpendic- 
ular to the easy axis. They studied the instability of 
domain wall structures for large external fields and gave 
variational stability boundaries including demagnetizing 
effects of fluctuations. 

In this paper we consider the different situation of an 
external field applied along the easy axis without the 
limitation to large external fields. The external field re- 
moves the degeneracy between the two anisotropy min- 
ima and consequently only pairs of Bloch walls can exist 
as static solutions. The basic topological excitations of 
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this system are thus twisted and untwisted pairs of tt- 
Bloch walls, r— I |— I 

Experimentstj and numerical simulationsll3 suggest 
that the annihilation of twisted domain-wall pairs in thin 
films requires much larger external fields than that of u: 
twisted domain-wall pairs. Furthermore, the observe' 
coercivity reduction in elongated particles at finite tem- 
peratures has no theoretical explanation. 

In this work, it is shown that both of these effects 
are related to the stability properties of twisted and 
untwisted domain wall pairs. The primary aim of the 
present paper is therefore a careful investigation of the 
fluctuations around these structures. We shall reveal the 
surprising fact that fluctuations around the twisted and 
untwisted domain wall pairs are described by the same 
set of operators. This puts the stability discussion of the 
untwisted and the twisted domain wall pair on an equal 
footing. It then follows immediately that the untwisted 
domain wall pair has exactly one unstable mode corre- 
sponding to an expansion or a shrinking of the structure. 
The untwisted domain wall pair is thus identifled as a 
"nucleus" of critical size in a first order phase transition 
and thus plays a cmpial role in thermally activated mag- 
netization reversaO in elongated particles. A detailed 
statistical mechanical theory of magnetiajation reversal 
will be presented in a forthcoming papeJiZI. 

Another immediate consequence of this relation is the 
instability of the twisted domain wall pair (or "27r-Bloch 
wall" ) for large extern£il fields as has been discovered by 
Magyari and ThomasEJ and independently in Ref. [ij. 
By a careful examination of the nonlocal demagnetizing 
fields which are not included in the model of a biaxial 
ferromagnet, it is shown that this effect should be ob- 
servable in thin films. In particular, the minimal attain- 
able distance of two domain walls is shown to decrease 
with increasing hard-axis ansiotropy. It is emphasized 
that this effect is beyond the otherwise highly successful 
descrigtion of domain walls within Slonczewski's effective 
modefl. The present results are also crucial, for the cur- 
rent design of vertical Bloch line memorieg23 whose read 
operations relytJ on a distinction between domain wall 
pairs with different relative sense of twist. 

The work is organized as follows. In section II we 
present the model and discuss its role as an effective 
model which describes planar structures in a 3D model 
including demagnetizing effects. In section III untwisted 
and twisted domain wall pairs are presented and their 
energy is evaluated. It is shown that both structures can 
be viewed as a coherent superposition of two 7r-Bloch 
walls. In section IV the operators governing the fluctua- 
tions around the 27r-Bloch wall and nucleus are derived. 
In section V these results are applied to discuss the in- 
stabilities of these structures. In section VI, bound state 
energies and scattering phase shifts of the fluctuation op- 
erators are discussed analytically and numerically in view 
of a calculation of nucleation rates of domain wall pairs. 
The discussion of scattering phase shifts provides a lu- 
cid example of the widely unknown version of Levinson's 



theorem in ID: Scattering phase shifts do not converge 
uniformly to those of the operators that are obtained in 
the limit of small and large external fields. In section VII 
we show that the present model can account for several 
different experimental configurations in thin films and we 
shall show that the nonlocal infiuence of demagnetizing 
fields on twisted domain wall pairs can be neglected for 
sufficiently thin films. The present model is thus ade- 
quate for small domain wall distances in sufficiently thin 
films where nonlocal demagnetizing effects are shown to 
be negligible. 

It is not necessary that the reader follows all details 
of the present paper. Those who are interested in exper- 
imental implications may skip the more formal sections 
IV and VI and directly proceed to section VII. 



II. THE MODEL 

In this we work we consider effectively one-dimensional 
magnetization configurations described by the following 
energy per unit area 



£ = 



El 



Mi 



(2.1) 



where M = M(z), 9^ = d/dz and A/q = |M| is the con- 
stan t magnitude of the magnetization. The first term in 
( ^.l|) is the classical counterpart of exchange energy and 
A is an exchange constant. The second term describes a 
hard-axis anisotropy characterized by the anisotropy con- 
stant Kh > thus rendering the xy-plane an easy-plane. 
The rotational invariance in this easy-plane is broken by 
an additional easy-axis anisotropy with anisotro py c on- 
stant Ke > 0. The last term in the integrand of (2.1) is 
the Zeeman term which is due to an external field Hoxt 
pointing along the easy- axis. 

Apart from the description of the (sublattic^nlspin con- 
figuration in ID (anti-) ferromagnetic affi|temsEH3, the en- 
ergy (2T) has found wide applicationstra in the descrip- 
tion of planar domain walls and their mobilities in bulk 
ferromagnets. As will be discussed in section VII, it is 
also adequate for the description of domain wall pairs in 
thin films. Due to the absence of discussion in the recent 
l iter ature, it seems convenient to review how the energy 
( ^.l[ ) may be derived from the energy of arbitrary 3D- 
magnetization configurations M = ]V[(r) in a volume V 
with inclusion of demagnetizing effects: 



E 



{VMyf + (VAf, 



(2.2) 



In contrast to (2.1), the first term in the integrand is 
the exchange term in three dimensions while the second 
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and third terms describe crystalhne easy- and hard-axis 
anisotropies of strengths isTe.cryst, ^ft,,cryst > 0, respec- 
tively. The fourth term is the demagnetizing energy 
with the demagnetizing field H„i obeying the magne- 
tostatic Maxwell equations V x = 0, V • B = 0, 
(B — H„i + 47rM). They can be rewritten in the form 
of a Poisson equation V^<j>m = 47rV • M with the mag- 
netostatic potential defined via Hm — — V$m- The 
Poisson equation is integrated in a standard way, and 
after splitting volume and surface terms we obtain 



H„(r) 



d^r' p,„(r') 



r — r' 



/|3 



+ / dS' a„,(r') 



r — r 



dV 



(2.3) 



where crm(r) = M(r) •n(r) is the magnetic surface charge 
(n is the normal of the surface dV) and pm(r) = V- M(r ) 
is the magnetic volume charge. Inserting (2_^) into (2.2) 
one recognizes that the evaluation of a mag netization 
configuration M(r) by minimization of (2_^) for given 
boundary conditions is in general a hopeless task. 

However, experiments reveal that the magnetization 
distribution in the vicinity of a domain wall in the bulk 
of a sample is a locally planar structure. This suggests 
the existence of an effective energy density which is of 
the form (|]^). In fact restricting ourselves to planar 
structures M = M.{z) and neglecting magnetic surface 



chargesEil in (^j) we obtain for an infinite sample a de- 
magnetizing field of the form 

H„(z) - -4^Af,(z)e„ (2.4) 

where is the unit vector in 2;-direction. For the deriva- 
tion of (2.4) we have also a ssum ed that Mz{ioo) = 0. 



After insertion of (2^) into (2^), the demagnetizing en- 
ergy takes the form of a hard-axis anisotropy along the 
z-direction. The underlying physical picture is simple: A 
planar arrangement of parallel dipoles has higher energy 
when the dipoles stick out of the plane than if they are 
in the plane. The form of the demagnetizing field ( p.4[ ) 
used to analyze wall motion experiments in garnet film; 
For stru ctur es of planar symmetry we thus may reduce 
(O) to (O) provided that 



K, = K, 



/i,cryst 



2ttM^ > 0, 



— -^e,cryst ^ 0- 



(2.5) 
(2.6) 



This holds for a configuration as e.g. shown in Fig. |^ a). 
For other sample geometries and anisotropy configura- 
tions, we can similarly express the effective anisotropy 
constants Ke, Kh in (2.1) by shape and crystalline 
anisotropies. 

To incorporate the constraint = Mq — const, in 



Eq. ( 2 . 1 ) , we use spherical coordinates defined by M /Mq 
= (sin cos 0, sin sin (/i, cos 9) . Further it is convenient 
to introduce dimensionless quantities by taking the scales 
of length and energy per area as 



M = [y] = [^] 




[£] = 2v/Ii^r. (2.7) 



Consequently, the units of the magnetic field are given 
by [H] = ^/W^. The length ^ AjK^ is the width of the 
static TT-Bloch wall, and 2\fAKl is half the energy per 
unit area of the static 7r-Bloch wall. With these defini- 



tions, the energy (2.1) becomes 



£ = j dz|-[(9,0)2+sin2 0(9,0)2] 
--[sin^ 6*0082 (/)-!] -I- -^0082 6* -/i sin 6*008 (2.8) 

where = 9{z) and = (f'i^)- The normalization is 
chosen such that the uniform states 9 = tt/2 and </> = or 
(j) — TT h ave zero energy in the absence of an external field. 
In ( ^.8[ ) we have introduced the dimensionless anisotropy 
ratio Q > with 



^ Kh' 



(2.9) 



describing the ratio of easy- and hard-axis anisotropy in 
the effective model (2.1). Note that this is a slight ex- 
tension of the common definition where Kh — 27rMg . In 
( ^.8[ ) we have also used the reduced external field h which 
is related to the external field i?oxt in laboratory units 

by 



2K, 



> 0. 



(2.10) 



At first sight, the choice of the coordinate frame in 
( p.lD and (2.8) might be surprising since the polar angle 
is not measured relative to the external field. The ad- 
vantage of such an orientation is that the linearization in 
the angles 9 and (j) around structures confined to the xy- 
plane is equivalent to a lincarizatio p i n .a cartesian frame 
that is rotated along this structuretiHl3 but is simpler in 
practice. Measuring 9 from the external field would not 
allow linearization in the azimuthal angle (j) to describe 
spin wave excitations of a uniform state parallel to the 
external field. 



III. DOMAIN WALL STRUCTURES 

In the following we s hall focus on static easy-plane 
structures of the model (2.8). It is shown that the only 
solitary easy-plane structures are twisted and untwisted 
pairs of 7r-Bloch walls. Simple representations are pre- 
sented that rel ate these solutions to each other. 

Inspecting ( ^.8|) , we recognize that the hard-axis 
anisotropy is minimized for 9 ~ 7r/2. The correspond- 
ing static structures then identically satisfy the Euler- 
Lagrange equation 5E/59 = while the Euler-Lagrange 
equation in (j) reads 



3 



sm cos ( 



h sin (h — 0. 



(3.1) 



Upon integration with d<j}/dz we obtain the first integral 
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with 



1 

^Kdz 



F(0) 



/i cos ( 



(3.2) 



(3.3) 



Eq. ( |3.2| ) has the form of an energy conservation for a 
fictitious particle moving in the one dimensional poten- 
tial V{4)). By this analogy, we can gain an overviewcJ of 
all static easy-plane structures. Note that the potential 
V{4)) is the negative of anisotropy and external field con- 
tributions to £ for = 7r/2 up to an irrelevant constant. 

Solitary solutions are now obtained as trajectories of 
the fictitious particle starting from a local maximum of 
V{(j)). Due to the "energy conservation" ( |3.2[ ) it will ei- 
ther creep into a different maximum of the same height 
or, if it started from a lower maximum, it will bounce 
back into the same state. For ft, ^ the degeneracy be- 
tween the maxima of V (i.e., minima of £) is lifted and 
two distinct trajectories emerge. One trajectory connects 
a global maximum oiV &i (j) — Q with an adjacent one at 
(j) = ±2^. This trajectory corresponds to a twisted pair 
of TT-Bloch walls. The second possible trajectory repre- 
sents a localized excursion from the lower maximum of 
V at (j) = ±7r which corresponds to an untwisted pair of 
TT-Bloch .walls. For other values of C in (3.2), periodic 
solutiongl3 occur which may be regarded as generaliza- 
tions of the above solutions to finite sample lengths. 

Thus we have gained an overview over all possible so- 
lutions without having solved the differential equation 
(3.2) in detail. This analogue should also prove useful 
for different models with other forms of the anisotropy 
and different orientations of the external field. 

Apart from the trivial symmetry arising from the 
repr esen tation of M in terms of spherical coordinates, 
Eq. ( |3.l|) is invariant under the following symmetry oper- 
ations 



and 



r : M -M, 

h —h. 



(3.4) 



(3.5) 



7?-x(7r) corresponds to an (internal) rotation of the mag- 
netization by an angle tt around the x-axis (0 = 7r/2), 
whereas T represents a time inversion. Therefore, all 
specific solutions quoted below have equivalents arising 
through the action of T?,^;, T and T o T^^. For a given 
direction of the external field there are thus exactly two 
equivalent structures related to each other by the action 



of 72,2: (tt). To classify the solutions it is also convenient 
to introduce the twist 



(3.6) 



Single TT-domain walls belong to \q\ = 1/2, whereas 
twisted domain wall pairs have \q\ = 1 and untwisted 
pairs have q = 0. Note that TZxiir) changes the sign 
of the twist whereas T leaves the twist invariant but 

reverses the magnetization at infinity. 

0, (E|) with (|3., 



For a vanishing external field, h = 0, (3.2) with ( |3.3| ) 
may easily be integrated with the boundary conditions 
920(±oo) — 0, 0(— oo) = TT, and 0(oo) = to yield the 
TT-Bloch wall 



bx{z) = 2arctane 
Ok = 7r/2. 



(3.7) 



The configuration (3/7) is shown in Fig. 0. The Bloch 
wall represents a smooth transition region between the 
two degenerate uniform states of minimal anisotropy en- 
ergy while the magnetization always lies in the easy 
plane. In (3.7) an integration constant describing the 
arbitrary wall position has been fixed such that the tt- 
Bloch wall is centered around the origin. However, as 
we shall see in Sec. IV, this degeneracy with respect to 
translations will lead to a (Goldstone-) mode of zero en- 
ergy in the excitation spectrum. The finite domain wall 
width arises through the balance of exchange energy and 
uniaxial anisotropy, the former tending to enlarge the 
transition region, the latter tending to narrow the Bloch 
wall. 

for ft = we obtain for the 



Inserting (3.7) into (2 



energy per unit area of the Tr-Bloch wall 



dz 



dz 



= 2, 



(3.8) 



where in the first step we have made use of the fact that 
4>K obeys the "energy conservation" ( [3.2| ) with C = 1/2 
and ft = 0. 

For ft ^ 0, the degeneracy between the two anisotropy 
minima at {9,(1)) = (Tr/2,Tr) and (Tr/2,0) is lifted. Con- 
sequently, single Bloch walls cannot exist any more. In- 
stead two different types of Bloch wall pairs arise which 
are discussed in the next two subsections. 



A. Untwisted domain wall pairs 

For < ft < 1, the boundary conditions ^^^(ioo) = 0, 
(f)(±oo) — TT, imply that C = lZ2=.ft. The integration of 
(|J) then yields the "nucleus"EaS 



<j)s{z) = 2 arctan 
9s = V2- 



/ coshz/^s 
\ sinhi?s 



(3.9) 
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As we shall see in Sec. V, the configuration (3.9) rep- 
resents a saddle point of the energy since it is unstable 
for all values < /i < 1 of the external field. Since it 
has exactly one unstable mode, it represents a critical 
nucleus f or m agnetization reversal. The integration con- 
stant in (3^) is chosen such that the symmetry center 
is located at z = . Note, however, that the continuous 
degeneracy of ( |3.9| ) with respect to translations will give 
rise to a zero energy (Goldstone) mode in the fluctuation 
spectrum, quite analogous to the case of the 7r-Bloch wall 
above. In ( |3.9| ) we have introduced the "radius" Rs of the 
untwisted domain wall pair. Rs is related to the external 
field h and the width (5, as follows 



h = sech Rs 



coth Rs 



(3.10) 



The nucleus may also be written as a superposition of two 
untwisted 7r-Bloch walls (|3.7|) centered at z/Ss = ±i?sj 



0,(z) = (j)K ( -J + Rsj + (l)K l^Y+Rs] ■ (3.11) 

Note that this relation is exact for all < i?^ < oo. Eqs 

:9i 



( |3.9[ ), (3.11) thus describe a domain with magnetization 
oriented parallel to the external field which is delimited 
by a pair of untwisted 7r-Bloch walls (cf. Fig. ||). The ex- 
istence of this structure is due to the balance of exchange 
and Zeeman energy. The exchange energy tends to at- 
tract the untwisted domain walls, whereas the Zeeman 
energy pulls them apart since it fa vor s the intermediate 
domain. As is illustrated by ( 3.10| ), ( ^.11 ), and Fig. ||, 
the domain wall separation tends to infinity for h —^ 
whereas for — > 1 the two oppositely twisted domain 
walls almost annihilate each other and the nucleus de- 
generates to an infinitesimal deviation from the uniform 
"down" state (0,6*) = (7r,7r/2). 
Using the parametrization ( 3.10| ), the energy per unit 



area (2.8) of the nucleus relative to the "down" state 



takes the simple form 



£s = £[(l)s,ds]~£[^ = 7r,9= -] = 
= 4 tanh Rs - 4i?s sech^ Rs. 



dz 



dz 



(3.12) 



In the first step we have used the first integr al (|3.2| ) and 
the integration is most easily performed with ( |3.11| ). The 



first term on the r.h.s. in ( 3.12| ) describes the deformation 
energy of the nucleus compared to the uniform state in 
the absence of an external field. The second term is the 
Zeeman energy —A4sh. The magnetic moment per unit 
area relative to the down state is thus given by 



Ms = iRs 



(3.13) 



The deformation energy vanishes for i?^ — > refiecting 
the fact that untwisted pairs of domain walls are attrac- 
tive. For Rs ^ oo, the energy converges to that of two 
independent 7r-Bloch walls. 



From (|1|) and (|1^) we may immediately derive the 
(formal) susceptibility 



dMs 



(3.14) 



This susceptibility has only formal character, since as we 
shall see below, the nucleus is unstable for all values of 
the external field < h < 1. 



B. Twisted domain wall pairs 

For < /i < oo, and for the boundary conditions 
d,(t){±oo) = 0, 0(±oo) = 0, we have C = 1/2 + h. 
Eq. ^luljiiay then be integrated to yield the 27r-Bloch 



4>biz) — 2 arctan 
Ob = 7r/2. 



/ coshi?{, 
\ sinh z/Sb 



(3.15) 



The integration constant has been chosen such that the 
symmetry center is located at z = but as in the case of 
the nucleus, the translational degeneracy will give rise to 
a zero energy (Goldstone) mode in the fluctuation spec- 
trum. The "radius" Rb of the twisted domain wall pair 
is related to the external field h and the characteristic 
width Sb as follows 



h — csch^ Rb 



Sh = tanh Rb 



(3.16) 



The 27r-Bloch wall ( 3.15| ) may also be written as a su- 
perposition of two twisted 7r-Bloch walls ( ^Tt] ) located at 
z/Sb = ±Rb, 



hbiz) = (t)K [ — Rb 



OK 



Rb 



(3.17) 



This relation is valid for all values of Rb- Eqns ( [3.15 ), 
( |3.17 ) describe a pair of 7r-Bloch walls located at z/S = 
±Rb with equal relative sense of twist, enclosing a domain 
of reversed magnetization (cf. Fig. ^). This structure is 
stabilized by the balance of Zeeman and exchange en- 
ergy. The Zeeman energy tends to enlarge the domains 
oriented parallel to the external field, whereas the ex- 
change energy pulls the twisted domain walls apart. As 
illustrated by Fig. |^ b), the 27r-Bloch wall decays for 
h ^ into two individual 7r-Bloch walls with increas- 
ing separation, whereas for /i — s- oo (Fig. ^ a)), the two 
TT-Bloch walls are squeezed and the transition region be- 
comes infinitesimally small. 

The energy per area of the 27r-Bloch wall is given by 

£b^£[(job,e=^]-£[cf> = 0,d=^]= r dz 



4 coth Rb + ARb csch^ Rb 



(3.18) 
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where (3.2) and ( |3.17 ) have been used. The first and 
second term on the r.h.s. in ( 3.18 ) describe the deforma- 
tion energy of the 27r-Bloch wall relative to the uniform 
"up"-state (/) = in the absence of an external field, and 
the Zeeman energy, respectively. The magnetic moment 
per unit area relative to the up state is thus given by 



Mt 



(3.19) 



Note that the deformation energy in ( 3.18 ) diverges for 
-Rfc ^ (i.e., h — > oo), i.e. a compression of the 2tt- 
Bloch wall to zero width is connected with an infinite 
increase in exchange energy. For i?^ — s- oo (i.e. h ^ 0) 
the deformation energy tends to that of two single tt- 
Bloch walls and the Zeeman energy becomes zero. With 
18) and (3.19) we obtain the susceptibility 



(I 



Xh 



dMb 



hy/T+h 



(3.20) 



For large external fields this susceptibility has a only for- 
mal meaning, since the 27r-Bloch wall can become unsta- 
ble for h^Q^^ /3 as we shall see in Sec. V. 



IV. FLUCTUATIONS 

To investigate the stability of the structures presented 
in the last section, we perform an expansion around a 
given easy-plane configuration {(f>Q{z), 9 = 7r/2) as follows 



4>{z) = (t>o{z) + <^(z), 
6{z)='k/2-p{z), 



(4.1) 



where |p| 1- First, we shall review the fluctuations 
around the 7r-Bloch wall because of their close relation 
to the fluctuations of the 27r-Bloch wall and the nucleus. 

for /i = we 



Inserting (4.1) with 0o = 'Pk into (^ 

obtain up to second order in (p{z) and p{z): 

1 1"^ 1 1'°° 

gi2)=£K^ ^^^^K^^ dzpin'' + Q-^)p, 

^ J —oo ^ J -oo 

(4.2) 

where is the Bloch wall energy (|3.8|) . No first order 
term in the fluctuations is present in ( [4.2| ) since (^k obeys 
the Euler-Lagrange equations ( p.l| ) with h — The 
operator H.^ is defined as 



K 



dz^ 



1 — 2 sech z. 



The potential appearing in (O) belongs to the family of 
reflectionless potentials which are of the form —m{m + 
1) sech^ z (m an integer), and which ar e di scussed in the 



appendix. The eigenvalue problem of (O) is 



iy = 0,k. 



(4.4) 



There is one bound state with zero energy 
1 



V2 



sech 2;, Eq 



K 



0, 



(4.5) 



and there are running (spin-wave) states 
1 



^27r(l + fc2) 
l + fc2. 



[—ik + tanhz] e 



ikz 



(4.6) 



The easy-axis an isot ro py le ads to t he g ap 1 in the spin- 
wave spectrum (4.5), (4.6) while (4.2) shows that the 
hard- axis anisotropy gives rise to the "mass" of the 
p-fluctuations. Since is nodeless and thus represents 
the ground state of Ti-K, all eigenvalues of TIk + Q^^ 
are positive. Therefore all fluctuations around a 7r-Bloch 
wall have positive energy except for the zero energy mode 
{(p,p) = {xo{z),0). This mode corresponds to a rigid 
t ran slation of the Bloch wall: Taking the derivative of 
( |3.2[ ) for h — we obtain Ti^^ d(f>K /dz = and therefore 
Xo oc d(j)K/dz. We conclude that in the absence of an 
external field the static kink is stable with respect to pla- 
nar distortions except for rigid translations which invohie 
zero energy. This result was first obtained by Wintertil. 
We now proceed with a discussion of fluctuations of the 
nucleus and the 27r-Bloch wall. 



A. Nucleus 



Inserting (4.1) with (f>o — (t>s into (2.8) and evaluating 
£[<l),0] — £[(f> = TT, 6* = ^] to 2"'' order in ip and p we 
obtain 

£(2) = + dz + - / dz pH'^Pp, (4.7) 

^ J-00 ^ J -00 



where Es is given by ( 3.12| ). The first order term in 
the fluctuations is absent since (f)s satisfies the Euler- 
Lagrange equations ( |3.l| ). The operators TL'^^ and 7i*^ 
are defined as 



= -—^ + 2 cos^ (j)s + sec\\ Rs cos (j^s-l, (4.8) 
dz'^ 



dz^ 



2cos'^ 



3 sech i?, cos ( 



2sech^ R,-l + Q-\ 



(4.3) with 



COS( 



sinh^ Rs — cosh^ (z/^s) 
smh^ Rs + cosh'^ i^/Ss) 



(4.9) 



(4.10) 



This form of the fluctuation operators is rather involved. 
Since the nucleus can be represented as a superposition 
( |3.11 ) of untwisted 7r-Bloch walls, we ex pect these oper- 
ators to contain potentials of the form (4.3) for each of 
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the constituents of the nucleus. In fact (4.8) and (|4 
allow for the much simpler representation 



with £b given by (|1|). The operators Ti.'"^ and are 
defined as 



■V- ( j,R, 



(4.11) W^*^ = -^ + 2cOs2(/)6 + CSch^i?6COS06 



(4.18) 



Q-\ (4.12) 



where the potentials V± are given by 

V±(C, i?) = 1 - 2 sech^(C + R)-2 sech2(C - R) 

±2sech(C + i?)sech(C-i?). (4.13) 



The second and third term on the r.h.s of ( [1.13| ) are the 
potentials ( |4.3| ) of two noninteracting 7r-Bloch walls lo- 
cated at z/db — -izRh- The last term, which vanishes 
for Rs ^ oo, describes the interaction of the two tt- 
Bloch walls and is thus sensi tive t o their relative sense of 
twist. The constant Q^^ in ( 4.12 ) is due to the hard-axis 
anisotropy and leads, in analogy to the 7r-Bloch wall, to 
a finite mass of fluctuations out of the easy plane. The 
corresponding eigenvalue problems are 



H^^xr (2, Rs) - (^s)xr (2, Rs) 
H'Px7iz,Rs) = EtP{Rs)x7i^,Rs)- 



(4.14) 
(4.15) 



where the index denotes bound states and scattering 
states. An analytical solution of these eigenvalue prob- 
lems seems only possible in the limiting cases i?^ — > 
and Rs oo. However, one bound state of H""^ , the 
zero energy state, can be derived immediately by tak- 
ing advantage of the continuous degeneracy of (g) with 
respect to translations. Taki ng t he z-derivative of (3.1) 
at = we obtain with (18) TC^f dcj) s / dz = 0, and 
therefore 



Sip 

Xi 



sech(— + Rs) — sech(- 



El 



Sip 



0. 



Rs) 
(4.16) 



The antisymmetry of the zero mode Xi''' with respect 
to z is a consequence of the opposite relative sense of 
twist of the two 7r-Bloch walls in (3.11). The remaining 
bound state energies and the scattering phase shifts will 
be investigated analytically and numerically in the next 
section. 



B. 2n Bloch wall 



Inserting (4T) with 0o ~ (pt into ( |2.8| ) we obtain for 
£[(t), e] -£[(l) = 0,e= f ] to 2"^ order in f and p 

£f ^ ^£b + ^ dzipn'"^ip + - dzp H'Pp, 

^ J — oo ^ J —oo 



(4.17) 



2 cos^ 



3 csch i?f, cos ( 



2 csch^ Rh-\ 



with (|t|) 



sinh^ (^/^b) ^ cosh^ i?b 
sinh^ {zjbb) + cosh^ Rh 



(4.19) 



(4.20) 



The operator (4.18) is identical to that describing the 
fluctuations around a kink in the do uble sine-G ordon 
modelEJ. In analogy to the nucleus, (4.18) and (4.1£) 
allow for a much simpler representation 



d^ 

dz^ 



-,Rb 



(4.21) 



H"" = (4.22) 



where the p otent ials V± are given by ( 4.13 ). The rep- 
resentation (4.21) has also been obtained by Sodano et 
aZ.Ej in the discussion of kinks in the double sine-Gordon 
mod el. It is instructive to compare (4.21) and ([4. 221) with 
(O): The 2"'' and 3''^ term on the r.h.s 



of 



are 



the potentials of the noninteracting domain walls located 
at z/5b — ±Rb- The last term describes the interaction 
of the two domai n wa lls and vanishes for i?f, — > cxd. The 
constant Q^^ in (4.22) is due to the hard-axis anisotropy 



and leads to a finite mass of out of easy-plane fluctua- 
tions. 

We write the eigenvalue problem of ( 4.21 ) and ( 4.22| ) 
in the following form: 

n'^^xtfiz, Rb) = Et^iRt)x'/iz, Rb), (4.23) 
n'^x'iriz, Rb) = EtP{R,)x''J'{z, Rb), (4.24) 

The index v denotes bound states and scattering states. 
Again, an analytic solution of these eigenvalue problems 
seems only possible in the limiting cases i?s — > and 
Rs — > oo. In analogy to the nucleus, one bound state of 
Ti.'^f can be derived immediately. Taking the z-derivative 



of ( 3T ) at (/) = we obtain T-L""^ d4> s / dz = and therefore 



dz 



Xo^ « ^ = \ sech(^ + Rb) + sech(f - R,) 



z 



e: 



(4.25) 



The symmetry of Xo^ with respect to z re flects the equal 
sense of twist of the two domain walls in ( |3.15| ) . 
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V. INSTABILITIES 



<o, 



(5.9) 



We are now in a position to st ate o ne of th e centr al re - 
sults of th is paper. Comparing ( ^.2l| ), ( f4.22| ) with ( |4.12 
and (4.11) we infer the remarkable connection 



Ob 



(5.1) 
(5.2) 



Here, for clarity, the notation Ti^'^{z,R) = —cP/dz^ + 
5~'^V-{j-,R), 5s = cothi?, has been used and analo- 
gously for the remaining operators. Eqns (5J^), (5.2) 
show that the fluctuations around the 27r-Bloch wall and 
around the nucleus are governed up to rescaling by the 
same set of operators. Consequently the eigenvalues are 
related by 



Et^R) 



Elf{R)+Q-\ 



and the eigenfunctions obey 

x7{z,R)^xlnj-^z,R), 

xl'',{z.R) = x7{^z,R). 

Ob 



(5.3) 
(5.4) 

(5.5) 
(5.6) 



where for bound states v = v' and for scattering states 
v = k, v' — {Ss/6b)k. The continuum eigenvalues are 



defined as El"^ 



v-2 



El' 



Q- 



- E-'T' for i 



In (5.1)-(5.6) we have used 



^b /5s — tanh^ R. 



(5.7) 



The relations (j^), ( |5!^ ) together with ( [1.16| ), ( [1.25| ) now 
allow us to discuss instabilities of the nucleus and the 27r- 
Bloch wall in a simple and straightforward way. 

The function x7 given in (4.25) is symmetric and 
nodeless, and hence it represents the ground state of Ti'^'^ 
with zero energy. Except for this state, Ti'"^ has a strictly 
positive spectrum and so has Ti"^, i.e. 



El^{R,)>Q, Ef,P{Rs)>Q, 



(5. 



for all I' and < R^, Rs < oo. It thus follows that i) the 
27r-Bloch wall is stable with respect to easy-plane fluctu- 
ations (neutrally stable with respect to the zero mode), 
and ii) that the nucleus is stable with respect to out of 
easy-plane fluctuations. 

On the other hand, the function Xi^ antisymmetric 
with one node and thus represents the first excited state 
of TC'^'^ . Since it has zero energy, there is exactly one 
nodeless, symmetric bound state of negative energy, i.e. 



for all < i?s < oo. The inequality (5J) is the origin of 
the following instabilities: 

nucleus: Fluctuations in (p direction exhibit exactly 
one mode of negative energy E^''' . Since H**^ is positive 
(cf. ( |5.8| )), we conclude that there is exactly one unstable 
mode of the nucleus for all values of Rs- Since (jjs is 
untwisted (i.e. q{(t>s) — 0), the instability in tp provides 
an example of a topologically induced instability. 

2TT-Bloch wall: The 27r-Bloc h wa ll is stable with respect 
to ip fluctuations because of ( |3.8[ ). Since q{4>b) = 1, this 
stability is of topological origin. However, an instability 
against out of easy-plane distortions occurs if 



E'o'\Rb) 



Q-^ -coth'^ Rb\E7{Rb)\ < 0. (5.10) 



Rb- 



wh ere R ^ is related to the exter nal field as ft, = csch 
In dslol ) we have made use of ( ^.4[ ). 

Anticipating results of the next section for the asymp- 
totic behaviour of the eigenvalues, we obtain the follow- 
ing asymptotic behaviour for this instability condition 



Q-^<2h, /i<l, 
Q-^<3h, /i>l. 



(5.11) 
(5.12) 



The numerically evaluated instability condition ( 5.10| ) is 
shown in F ig, m together with its asymptotic behavior 
( ^.llD and ( |5.12| ) . The instabihty of the 2tt Bloch waU|iS| 



in accordance with the result of Magyari and Thomasta 
who gave also an improved analytical estimate of the in- 
stability range for large h, however they did not discuss 
the nucleus and the relation of its fluctuations to the 2tt- 
Bloch wall. 

Since we have shown in this section that the eigenfunc- 
tions and eigenvalues of Ti^P and Ti^'^ can be expressed 
by those of Tf^ and H""^, we may restrict ourselves to a 
discussion of the latter operators in the following. 



VI. DISCUSSION OF AND H"'' 

In this section we evaluate the eigenfunctions of W"^, 
Ti.'^P numerically and provide analytical results in the lim- 
its of large and small Rs- We first discuss bound state 
energies which are related to the stability properties of 
the 27r-Bloch wall and the nucleus. In view of statistical 
mechanical approximations, the scattering phase shifts 
of the continuum eigenfunctions are discussed. Further- 
more it is shown that the appearance of zero energy res- 
onances in the spectrum require a subtle analysis of the 
applicability of analytical approximations. 



A. Bound states 

In the limit of large and small Rs , the eigenvalue prob- 
lems of H^''^ and can be solved exactly: 
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For large Rg, the potentials V±{j-, Rg) decay into two 
independent wells of the form — 2^J^sech^(^ ± Rg) and 

we denote the corresponding operators by Ti.'^ . This limit 
of large R^ is sometimes also referred to as "thin-wall 
limif'Ea. The bound states of H'"^ and are then 
given by the symmetric and antisymmetric combinations 
of the bound states of the single wells. For i?s — > cxi we 
thus have Xo^ ~^ Xo^ a-nd xT Xi^i where 



Xo'^(z) cx sech(-^ + Rs) + scch(-^ - Rs), 

Os Os 

Z Z 

xl''{z) cx sech(— + Rs) - sech(- Rs). 



(6.1) 
(6.2) 



Note that the r.h.s. of (6.1) and ( |6.2| ) are the exact zero 
energy eigenfunctions of Ti.^^ — and TL"''' , respec- 
tively. Since for large Rs these operators differ by a term 
©(e"^^"), we obtain within first order perturbation the- 
ory: 



E^iRs) 

- s-' 

~ -8e- 
and 

EfiRs) 



{x7,n^^x7) 
(xr,xr) 

3 1 



cosh Rs 

2R^ 



2Rs - sinh 2Rs 



sinh^ Rs 2Rs + sinh 2Rs 



(6.3) 
(6.4) 



(xf,H-^xf) 
(xf,xf) 

1 



2Rs + sinh 2Rs 



sinh i?< 



Q 



cosh Rs '^Rs 

-1 



sinh 2_R, 



Q-\ (6.5) 
(6.6) 



where (u, v) denotes the standard scalar product J dzu*v. 

For small Rs, we have H""^ ^ Ti^"^ and H^p 
with 



H'P = - 



d^ 

dz^ 

dz^ 



S-^[l - esech^ 



Ss^[l - 2scch2 



(6.7) 
(6.8) 



Both potentials (|6.7j ), ( |6.8[) belong to the class of reflec- 
tionless potentials which are discussed in the appendix. 

The (unnormalizcd) bound states of H'^'^ and their en- 
ergies are given by 



xl'^iz) = sech^ 



^0 



'36: 



Xi^iz) = sech-^ tanh E^'^ = 

Os Os 



0, 



and the spin wave states read 



Xh'^iz) = ( 3 tanh^ -j SikSs tanh -r 1 — (kSs 

\ Os Os 

E7 = + k\ 



(6.9) 
(6.10) 



*2 \ ikz 



(6.11) 



The op erat or H^p is up to rescaling analogous to the 
operator (O) which describes the fluctuations around a 
single TT-domain wall. It has one bound state Xo'{z) = 
sech{z/ds) with energy E^p = Q^^, and spin wave states, 



(6.12) 



Xk{z) = y-ik6s + tanh - 
ElP ^ Q-^ + 5-^ + k\ 



In Eqs (3.7)-( p.l2| ) we have to put 5s ~ Rs in order to be 
consistent with the terms neglected in the derivation of 
WP and W'P. 

We are now in a position to verify the asymptotic be- 
havior of the in stab ility threshold of the 2 7r-Blo c h-wal l. 
Inserting (0), (|j|) into ( |5T0| ) we obtain (|]Tl|), ( p^ , 
respectively. 

For arbitrary Rs , the bound state energies of H^'^ and 
Ji^P have been evaluated numerically and the results 
are summarize d in Fig. Jq. The values of the asymp- 
totic formulas (3.3) and ( |6^5D are represented by dashed 
lines. Note that they are accurate for values as small 
as Rs — 1.5. The operator Ti*"^ has three bound states, 
the ground state of negative energy Eq'^ , the zero- mode 
with E'l'^ = and a weakly bound state whose energy is 
always within 1% of the continuum threshold according 
to numerical calculations. For applications such as the 
evaluation of nucleation rates we can therefore use 



(6.13) 



This bound state does not seem to be a numerical artifact 
since its existence also follows from the long-wavelength 
behaviour of the scattering phase shifts A^^^ as we shall 
see in the next section. The ground state wave function 
Xo^ can be considered as an internal "breathing" mode 
of the nucleus, corresponding to an expansion or shrink- 
ing, dependin g on the sign of Xg"^. Note, however, that 



according to (6T), a strict equality Xo*^ 



14is/dRs only 

holds in the limit Rs ^ oo. The operator T-^p always 
has two bound states. The ground state with constant 
energy EqP = > has its origin in the Goldstone 
mode of the 27r-Bloch wall while the excited state Xi'' 
of energy i?*^ is related to the "breathing" mode of the 
27r-Bloch wall. 

Comparing the previous analytical discussion with 
these numerical results we are left with a paradox. Ti^"^ 
and both have two bound states, whereas numerical 
calculations reveal the existence of three bound states of 
Ti^'^ . Similarly, 1-i^P has one bound state whereas li^ and 
Ti^P have two bound states. The resolution of this para- 
dox lies in the fact that each of the operators obtained in 
the limits i?s — > 0, cxd exhibits a zero energy bound state. 
Any increase in the potential strength thus leads to an 
additional bound state which is precisely the reason for 
the excess bound states of TL"''^ and WP . The two well 
approximation has the same number of bound states 
as TC^P but an additional zero energy resonance. As Rs 
becomes finite, the zero energy resonance of Ti'^ is shifted 
into the continuum. 
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B. Scattering Phase Shifts 

The knowledge of scattering phase shifts is of impor- 
tance for statistical mechanical applications. In par- 
ticular, the results of tiw present section will be used 
in a forthcoming articlet3 on nucleation of domain wall 
pairs. The scattering phase shifts A^^^ (^(e)) ^^'^ 
(even) eigenfunctions Xk(o)^ (Xk(e)) of the operators 'W", 
i = ip,p are defined as follows: 

X^^(e) {z ^ ±00, Rs) cx cos [kz ± A^^) (fc, i?,)/2] , (6.14) 
xl\^^iz^±^,R,)^sin\kz±All^{k,R,)/2\ , (6.15) 



where i = ip,p. It is sufhcient to restrict our considera- 
tions to the phase shifts of W"^, since according to (5.5) 
(Em we have 



A;i;(fc,i?) = A-(|fc,i?), 

A%ik,R)^A^/fk,R), 



(6.16) 
(6.17) 



where j — e,o. 

For large i?s, the potentials —26~^sech^(z/ds±Rs) act 
as independent scattering centers, each contributing a 
phase shift 2 arctan(l/fc(5s). Therefore we have 



A%k) = 4arctan 



1 

k6. 



(6.18) 



For small i?^, the continuum j^genfunctions (6.11) and 
( |6l2| ) of n^'P and H'p lead toEZi 



A'^^ik) = 2arctan 



3kSs 



{kSs^ 



A'Pik) = 2arctan— . 



(6.19) 
(6.20) 



Eqns. ( 6.18| )-( 6.20| ) do not distinguish between odd and 
even parity eigenfunctions. 

It is a surprising fact that some of the scattering phases 
= (p,p, j — e,o do not converge uniformly to 



A 



111 



(5.18-6.20) in the limits Rg —>■ and Rs 00, respec- 
tively. Numerical calculations show that (cf. Figs. 

A^^^)(fc^O,i?.) =3^, 



A^^)(fc^O,i?.)=^, 
A^^)(fc^0,i?,)=2^, 



3.21) 



for all Rs- Eqns ( |6.18 )-(6.20), however, deliver the rela- 
tions 



This discrepancy has the same roots as the paradox en- 
countered in the previous subsection, namely the exis- 
tence of zero energy resonances. This is elucidated bj4 
the widely unknown ID version of Levinson's theoremES 
which relates the long wavelength limit of the phase shifts 
with the number of bound states: 

The odd-parity wavefunctions behave like in a 3D scat- 
tering problem: 



A^:)(fc->0) = 2^iV, 



(6.23) 



where N^^-^ is the number of odd-parity bound states of 

W**. In the absence of zero energy resonances, the scat- 
tering phase shifts of even-parity wavefunctions obey a 
different relation 



1, 



A^:)(fc^0) = 2^(iV(*,)--), z = ^,p 



(6.24) 



where iVj*^-) is the number of even-parity bound states of 

Ti'*^ According to the previous subsection, we have for 
all values of i?» 



(o) 



1, 



(6.25) 
(6.26) 



This shows that the A; ^ behavior of the scattering 



phase shifts (6.21) is in complete agreement with the 
number of bound states of Ti^"^ and as evaluated 
in the previous subsection. 

If zero energy resonances are present and if the po- 
tential i^reflectionless, Eqs. (3.23), (6.24) have to be 
replacedE3 by the parity independent expression 



Al}^{k^O)=n(NI^)+N, 



'(e) 



(6.27) 



where i = (p,p and j — e,o. Eq. ( 6.27 ) relates ( |6.22 ) to 
the number of bound states of the reflectionless operators 
H'"^, n^P, and as given in (|6^-(|6.12D. Levinson's 
theorem thus relates the nonuniform convergence of the 
scattering phase shifts towards A'"'', A'^p, A'^ to the ap- 
pearance of zero energy resonances in H'^'^ , 7Y*p, and Ti.'^. 
How this subtlety affects statistical mechanical consider- 
ations, will be discussed in a forthcoming papeiO. 

The short wavelength behaviour of the scattering phase 
shifts can be described within Born's approximationcj. 
For an operator —d^/dz"^ -\- V{z) with a symmetric po- 
tential {V{z — ^ ±00) = 0), the phase shift is given by: 



tan - 



1 

'2k 



dzV{z)sin^{kz). (6.28) 



If fc ^ is much smaller than variations in V{z), we can use 
sin^( fcz) ~ 1/2 a nd af ter insertion of V±[z/5sj Rs) — 
(cf. ( |4.13| )) into ( |6.28| ) we obtain 



A"'^(fc- 


-0) 


= 2^, 




Ag)(fc) 


A'^(fc- 


-0) 


= 2^, 




2 


A'P{k- 


->0) 


= TT. 


(6.22) 





1 r„ Rs 

2 =F 

k5s sinh Rg cosh R, 



], for kSs > 1, 

(6.29) 
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where the upper sign refers to the i = p and the tan- 
function has been replaced by its argument. Finally it 
is int erest ing to note that for H^'p the Born approxima- 
tion ( 6.28 ) with sin^(fcz) ^1/2 coincides with the exact 
result (|6.2C|). 



VII. 2n BLOCK WALLS IN THIN FILMS 

The results of the previous sections are rigorous within 
the ID model of a biaxial ferromagnet which contains ex- 
change, local anisotropics and the coupling to an external 
field. While we have seen in Sec. II that local demag- 
netizing effects can be incorporated into the model by a 
redefinition of the anisotropy constants, one might ques- 
tion the applicability to thin films where nonlocal de- 
magnetizing effects are not a priori negligible. Since the 
nonlocal demagnetizing interaction decays algebraically 
while the exchange interaction between domain walls de- 
creases exponentially, we expect the exchange interaction 
to be dominant and thus our model to be adequate for 
small domain wall separations. 

Indeed, it is the purpose of this section to show that 
for sufficiently thin films and at external fields below 
the threshold ( 5.10 ), twisted domain wall pairs may be 
brought sufficiently close such that the exchange inter- 
action between the individual domain walls becomes im- 
portant and nonlocal demagnetizing effects become irrel- 
evant. In this case our model adequately describes the 
equilibrium separation of the walls. We shall use c.g.s 
units throughout this section. 

To be specific, we choose coordinate axes as in Figs. 
^ a), b) and consider a film of thickness D in x- 
dircction which extends infinitely in the y-direction 
and has length L in ^-direction. Further we assume 
the magnetization to be strictly one-dimensional, i.e. 
M = ]V[(z). The demagnetizing energy per area f„i = 
— (1/2D) J dx dz llm{x, z) ■ M(z) can then be cast-into 
a very convenient form due to Dietze and Thoma£3: 



.L/2 

£m= / dz2TTM'^{z) + 



1 fL/2 

- / dz dz' [M,(z)Af,(z') - M,{z)M,{z')] x 



D 
In ( 1 



L/2 J-L/2 

{z-z'Y 



Eq. ( 7.1 ) reduces to a simple hard-axis anisotropy energy 
in the following two limiting cases: 

For a film thickness D smaller than the characteristic 
length scale of M, i.e. smaller than the domain wall 
width, the integration over t he re lative coordinate in the 
second term on the r.h.s. in ( [7.l| ) can be performed, and 
Eq. (7.1) reduces to 



L/2 



L/2 



dz MI{z). 



(7.2) 



This has the form of a hard-axis anisotropy energy nor- 
mal to the film plane. 

In the opposite limit of large D ^ L. the second term 
on the r.h. s. of (7.1) tends to zercO and we recover 
the result (2.4). This form of the demagnetizing energy 
is used for the description of domain wall dynainics in 
moderately thin (~ 1/im) rare earth garnet filmsElEl. 

As the magnetostatic interaction has the form of a 
hard-axis anisotropy in these limits, the energy density 
( ^.8| ) can be used to describe three distinct experimental 
configurations. In addition to the bulk situation consid- 
ered so far (cf. Fig. ^ a), it describes configurations in 
thin films {D / A/ Kf. < 1) which are perpendicularly 
(Fig H b) or in-plane (Fig || c) magnetized provided that 
the coordinate axes are chosen appropriately. The results 
of the previous sections thus hold for all configurations 
shown in Fig. || in the limit of infinitesimally thin films. 
In the following it is discussed how nonlocal demagne- 
tizing fields, i.e. the nonlocal contribution in (7.1), will 
affect these results for a film of finite thickness. 



A. Perpendicularly Magnetized Films 

Consider a situation as in Fig. ^ b) which requires 
the crystalline easy-axis anisotropy to be larger than the 
demagnetizing energy, i.e. i^e.cryst > ^ttMq, and the 
easy-axis to be oriented perpendicularly to the film, a 
situation typically realized in bubble films or in Barium- 
Ferrite. To estimate the nonlocal demagnetizing inter- 
action for domain wall separations large compared to a 
domain wall width, we consider the configuration 



M(o)(r) = 



MqBx for L > \z\ > a 
'Moex for \z\ < a. 



(7.3) 



for 0<a:;<r',— oo<y<(X) and vanishing else- 
where. The nonlocal magnetostatic interaction is due 
to the nonuniform surface charge distribution caused by 
the reversed domain at \z\ < a (cf. Fig |^ b). Inserting 
(7.3) into ( [7.l| ) and performing the Uffiit L oo after 
evaluation of the integrals, we obtainEil 

f Pr'' = f„[M(o)]-£,„[Moe,] = 



AM^Da I 4 tan 



(7.1) 



- -H2alna + {- - a)ln(l-fa^) } , 
a a ' 



where a — 2a/ D is the width of the reversed domain with 
respect to the film thickness. Note that the r.h.s of (7.4) 
decreases with increasing a and thus favors an expansion 
of the reversed domain independent of the relative twist 
of the domain walls. For the twisted domain wall pair it 
thus competes with the repulsive exchange interaction. 

An external magnetic field i?oxt along the positive x- 
dircction (cf. Fig. 8 b) will counteract this magnetostatic 
repulsion. A relation for the corresponding eq uilibrium 
width is obtained by minimizing the energy ([7.4[) plus the 
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Zeeman energy of the intermediate domain, 2Mo-ffcxt-Da, 
with respect to a with the result 



Hext 2 

— rr = ~ tan 
AttMq tt 



1 hi 1 + — 

a TT \ a'' 



(7.5) 



Note that the relation ( |7.5D does not depend on the ab- 
solute size of the reversed domain but only on its relative 
size with respect to the film thickness. Since on the other 
hand the mutual exchange repulsion of twisted domain 
walls depends on their absolute distance, it can only be 
observed if the equilibrium width in Eq. ( |7.5| ) is small. 
However, since the external field must not exceed the 
instability threshold, this can only be achieved in suffi- 
ciently thin films. 

To investigate this effect quantitatively, w e ha ve to 
compare the equilibrium width 2a = aD of ( [7.5[) with 
the separation 2a = 2R b SbSQ o f the domain walls forming 



a 2tt- Bloch wafi (|3.16| ), ( |3.17| ) where So = y/A/Ke is the 
static Bloch-wall width. It is sufficient to look at small 
HoxtMo/ (2Ke) where this relation between equilibrium 
width and external field can be expressed as 



H, 



4-kMo 



2K, 
ttM§ 



exp 



2a 



(7.6) 



This topological (exchange) interaction between the 
domain walls thus decreases exponentially, while the 
magnetostatic interaction ( |7.5| ) decays algebraically, i.e. 
Hcxt 1 47rA/o = Dj (27ra), for large 2a/ D. Since the r.h.s of 
( [7.6D is proportional to the "quality factor" Ke/{2iiMQ) 
of the domain wall, the exchange repulsion will manifest 
itself at larger domain wall separations with increasing 
quality factors and decreasing film thickness. 

This is illustrated in Fig. |9| where the so lid l ines show 
the exchange dominated w all separation (7.6) whereas 
the dashed lines represent ( [7.5[ ). The experimentally re- 
quired field to reach a certain wall separation will follow 
the curve that has the maximal value of i/cxt consistent 
with the sample parameters. Note that the exchange re- 
pulsion between the domain walls has a drastic effect. 
E.g. for a film with D/5q = 0.1 and K^/{2ttM^) = 10, 
the exte rnal fields to esta blis h a distance of 75 predicted 
by ( [7.6D exceeds that of ( [7.5D by an order of magnitude, 
or vice versa, the wall separations differ by a factor of 3 
for i7oxt/(47rMo) = 0.01. This large discrepancy should 
be accessible to experimental verification. 

So far we did not discuss the case of untwisted domain 
wall pairs. Despite the fact th at t hey have been shown 
to be unstable in Sec. V, Eq. ( [7.5| ) shows that such do- 
main wall pairs can exist in thin films due to the balance 
of demagnetizing and Zeeman energy provided they are 
well separated. However, the ansatz (7.3) overestimates 
demagnetizing effects. As the wall separation decreases, 
the untwisted walls annihilate each other (cf. Fig. |^ 
a)) and the magnetostatic surface charges are drastically 
reduced compared to those of (7.3). This implies that 



pairs will not follow the dashed curve in Fig. || down to 
vanishing a but exhibit an instability at finite a. This 
instability will occur at fields that are much smaller than 
the instability threshold of the twisted domain wall pair. 
We do not consider this case further but conclude with 
the remark that a quantitative theory can be obtained by 
improving the ansatz (7.3) by replacing it by the nucleus 
solution {<j)s, 9s). 



B. In-Plane Magnetized Films 

Consider the configuration shown in Fig. ^ c) with 
a crystalline easy-axis in the film plane. All previous 
formulae hold also for this configuration provided that 
we redefine all coordinates appropriately, i.e. {x,y,z) 
{z,x,y). Since the magnetization configuration is exclu- 
sively in the film plane, there are no induced magneto- 
static surface charges in a film that extends infinitely 
in x-direction and there is no equilibrium domain for- 
mation in the infinite film geometry. However, we as- 
sume that domain walls exist that have been created 
e.g. by nucleation at a sample end and/or nonuniform 
external fields. If domain walls are present, a magneto- 
static interaction arises between the magnetostatic vol- 
ume charges of the domain walls. For suffi cient ly well 
separated twisted domain wall pairs, Eq. ( |3.17 ) takes 
the form My/Mo = sech((?/ + a)/So) — sech{{y — a)/6o). 
Inserting this into (7.1) (with redefined coordinate axes) 
we obtain for a much larger than a domain wall width 
So, 



Sir, - 2D{7TMof 



(7.7) 



Eq. ( [7.7[ ) is simply the interaction energy of two infinitely 
long strings of dipoles along y with dipole moment per 
unit length = D J dyMy{y) = DSot^Mq- The mag- 
netostatic interaction between the walls is thus repulsive 
for the twisted domain wall pair and competes with the 
exchange interaction between the individual walls. (For 
untwisted domain wall pairs the magnetostatic as well as 
the exchange interaction would be attractive.) An ex- 
ternal magnetic field in x-direction will counteract this 
repulsion. The equilibrium distance 2a between two do- 
main walls is obtained by minimization of demagnetizing 
plus Zeeman energy, £^ -|-4Afo77oxta, with the result 



Ho 



AttMo 



■K D 
2T0 



the experimentally observed separation of untwisted wall 



In Fig. ^ b), this is compared with the wall separation 
(7.6) which is predicted by our model. Demagnetizing 
effects are obviously weaker than in a perpendicularly 
magnetized film. For films of thickness D = Q.\5o and 
Ke/ {2ttM'^) = 10, the exchange interaction dominates 
the demagnetizing interaction already at a domain wall 
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separation of 2a = 13(5o which is surprisingly large con- 
siderin g th e exponential decrease of the exchange inter- 
action (^^). Note also that the external fields which are 
required to achieve a domain wall distance of 6(5o differ 
by a factor of 100. Finally, we note that untwisted wall 
pairs in in-plane magnetized films are never stable for the 
anisotropy configuration shown in Fig. |§ c). 



VIII. CONCLUSION 

In this paper we have discussed the stability of twisted 
and untwisted domain wall pairs within a ID model of 
a ferromagnet. The fluctuations around these structures 
have been shown to be described by the same operators. 
By means of exactly known eigenfunctions which are re- 
lated to the Goldstone modes of these structures it has 
been shown that untwisted domain wall pairs exhibit ex- 
actly one unstable mode while twisted domain wall pairs 
are subject to an instability at large external fields. 

Furthermore we have argued that this model ade- 
quately describes the separation of twisted domain walls 
in ultrathin films and thus the above instability should 
be observable. 

Although untwisted domain walls are unstable within 
the biaxial ferromagnet, they can exist in thin films at 
large separations due to the long range magnetostatic in- 
teraction. However, as a consequence of their topological 
instability, the corresponding collapse field will be much 
smaller thanjthat of twisted dcpaain wall pairs. There are 
experimental and numericaO hints for this behavior, 
but more systematic studies are required to allow for a 
quantitative comparison with the present theory. 

Another important aspect is the role of untwisted do- 
main wall pairs as nuclei for magnetization reversalriji 
elongated particles. As has been reported elsewherella, 
the existence of such nonuniform nuclei fact can lead to a 
substantial reduction of the coercivity compared to stan- 
dard theories of magnetization reversal. Further details 
of the statistical mechanical theory of rnagnetization re- 
versal are covered in the following papeiO. 
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APPENDIX 



In Q, Q, (I 
ators of the form 



we encountered Schrodinger oper- 



m(jn -f 1) 
cosh^ X 



(Al) 



with m an integer. In the following, we shall show how 
the corresponding eigenvalue problems may be solved ex- 
actly with the method of Ref. The continuum eigen- 
functions of TiS^^^ have the remarkable property that 
their reflection co effic ient is zero. We write the eigen- 
value problem of (Al) as follows 



(A2) 



The key point for the solution of the eigenvalue problem 
( [A2| ) is the observation that TY'^™-' may be factorized in 
two different ways: 



= Q'"'+'^Q^"+'^-(m + l)2, 



with 



ax 



; tanh c 



Operating on ( |A2| ) from the left with Q^"', Q 

we recognize that if V'i™'' eigenfunction of 7Y*^™' with 
eigenvalue A, then 



(A3) 
(A4) 

(A5) 

{m+l) 



, (m+l) ^(in+1) , (m) 



(A6) 
(A7) 



are eigenfunctions of H^"^^^\ TL^"^^^\ respectively, to 
the same eigenvalue A. 

The continuum eigenfunctions of 7i("') can thus be re- 
lated to those of H^"^~^\ Since 7^^°) represents the free 
problem, the continuum eigenfunctions of Ti^™) can be 

obtained by successive application of Q^"' onto plane 
wave solutions, i.e. 



Jm) 



(A8) 



and belong to the eige nvalu e A = fc^. For Ti'-'^^ 
-cf/dx'^ - 2sech\, Eq. yields 



\-ik + tanh xje*'''' 



(A9) 



For HP'^ = —(P/dx^ — Gsech^x we obtain the contimmm 
eigenfunctions 

V'f = QV'^e''^ (AlO) 
= [Stanh^a; - 3ifctanha; - 1 - /c^jgifex^ j-^-,^-,^-) 

To find the bound states with A < 0, we first remark 
that the normalization of bound state eigenfunctions with 
different m are related as 
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+ CX3 



+00 



= (A + m^) / dxi^i"'^^. (A12) 



Continuing this recursion towards lower values of m we 
recognize that the normalization of, say, V-'i' would be- 
come negative. This can only be avoided if the recursion 
( A12| ) stops, i.e. if the bound state eigenvalues are given 
by 



A/ 



-I' 



1 = 1,2,. 



(A13) 



According to ( A12 ) this implies that 

0. (A14) 
This differential equation can be integrated with (|A5|) 



(x) = sech X. 



(A15) 



For m > I, the unnormalized 1*'* bound state (counted 
from the continuum) can be obtained recursively with 
the help of (M) 



I (m) 

V'a =q+ 



Specifically, we obtain for m = 1 
= sechx, 

with energy Ai = —1, and for m = 



(2) 



■02^'' = sech^x, 
sechxtanhx. 



(A16) 

(A17) 

(A18) 
(A19) 



with energies A2 = —4 and Ai = —1. All operators (Al) 
share the property of having a zero energy resonance. 
This means that an infinitesimal increase in the poten- 
tial strength of (Al) leads to an addition al b ound state. 
Therefore, the occurence of the operators ( Al) as describ- 
ing fluctuations around nonlinear structures in some limit 
of the external field has to be handled with care, since 
their number of bound states in general differs from those 
of the operators they emerge from. 
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FIG. 1. a) The 7r-Bloch wall interpolates between differ- 
ent anisotropy minima; b) Fluctuations tp, p around a given 
structure with (j)s,ip and 6 = n/2 at a given space point z. 
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FIG. 2. The nucleus is shown for a) small fields {Rg — 3.5) 
and b) for fields close to the anisotropy field {Rs = 0.4) 



FIG. 3. The 27r-Bloch wall is shown for a) small fields 
{Rt = 3.5) and b) for large fields {Rt = 0.4) 



FIG. 4. Stability and instability regions of the 27r-Bloch 
wall as a function of the external field h and the demagne tizing 
field strength . The dotted and dashed lines refer to (5.11) 
and (5.12) respectively. 



FIG. 5. The rescaled bound state energies of H"'^ and W 
are shown as a function of Rs {Ss = coth_Rs). The shaded 
region indicates the continuum states. The horizontal lines 
El'^ and Eq" correspond to the zero modes of H""^ and n''^, 
respectively. i?2^ is a very weakly bound state just below the 
continuum threshold. The dashed lines indicate the approxi- 
mation formulas (6.3) and (^.5|). The bound state energies of 
H'"^ and H''^ may be obtained from (0),(Ei). 



FIG. 6. The odd and even parity scattering phases of Ti.'"'' 
are shown for different values of Rs {5s = coth_Rs). A**'(fc) 



and A''(fc) are given by (3.18) - (6.2C). Note that the conver- 
gence A^^^j(fc) ^ A^^ik) for Rs^O and A^^^ik) A=(fc) for 
-Rs ^ oo is nonuniform. 



FIG. 7. The odd and even parity scattering phases of Ti."^ 
are shown . The convergence Aj^j(A;) A''*'(fc) for Rs 

and A^^j(fc) A"{k) for Rs ^ cx) is nonuniform. 



FIG. 8. Magnetization co nfigurations in films that can be 
described by the model (2.8). Note that the anisotropy con- 
stant has to be chosen as indicated to incorporate the local 
part of the magnetostatic interaction. Note also that the ori- 
entation of the coordinate frame in a) as used in the text is 
different from b) and c). 



FIG. 9. Domain wall separations for a twisted domain wall 
pair in a) perpendicularly and b) in-plane magnetized films 
due t o ba lance between external field and exchange (solid line, 
Eq. ( |7.6| ) ) or betwee n ex ternal field and demagnetizing effects 
(broken line, a) Eq. (ItJ)), b) Eq. (Ei). 
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